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ABSTRACT 

Discrimination  between  underground  nuclear  explosions  and 
shallow  earthquakes  using  the  vertical  component  of  the  Rayleigh 
wave  is  achieved  using  several  techniques.  This  analysis  includes 
the  previous  developed  area  under  the  Rayleigh  wave  (ARZ) ,  the 
newly  applied  total  energy  (ERZ),  and  the  total  energy  transported 
across  a  unit  width  of  the  waveguide  by  the  Rayleigh  wave  (REF) . 
Results  for  several  explosions  and  earthquakes  of  varying  magni¬ 
tudes  are  presented.  Evaluation  cf  these  techniques  and  their 
applicability  in  an  automated  discrimination  program  is  discussed. 
An  attempt  is  made  to  incorporate  a  matched  filter  approach  for 
weak  signals. 


INTRODUCTION 


The  objective  of  this  study  is  to  examine  possible  discrim¬ 
ination  techniques  using  the  vertical  component  of  the  Rayleigh 
wave.  The  use  of  surface  waves  in  a  discrimination  process  was 
originally  discussed  by  Brune,  Espinosa,  and  Oliver  (1963)  with 
their  introduction  of  the  AR  diagnostic  (the  area  included  in 
the  envelope  of  the  surface  waves  of  a  seismogram) .  Subsequent 
investigation  has  shown  that  better  results  are  obtained  using 
only  the  vertical  component  of  the  Rayleigh  wave  in  the  AR  pro¬ 
cess.  With  this  tool  as  the  basic  discriminant,  and  with  the 
idea  of  obtaining  optimum  results  for  the  time  invested,  modi¬ 
fications  have  been  made  to  the  AR  process,  and  new  discrimi¬ 
nants  have  been  tried.  A  program  has  been  written  which  yields 
results  for  each  discriminant,  enabling  comparisons  to  be  made 
readily.  In  the  following  paragraphs,  a  theoretical  description 
is  included  for  each  tool,  the  results  for  each  discriminant  on 
actual  data  is  given,  and  the  program  which  performs  these 
calculations  is  described. 

Methods  of  Analysis 

ARZ :  This  was  originally  defined  as  the  area  within  the 
envelope  of  the  vertical  component  of  the  Rayleigh  wave  recorded 
on  a  standard  instrument  at  standard  magnification.  This  quantity 
has  no  physical  meaning.  Its  motivation  was  founded  on  the  fact 
that  it  was  easy  to  measure  manually  by  tracing  out  the  envelope 
and  then  using  a  planimeter  to  obtain  the  area. 

In  our  effort  to  automate  this  process  and  bring  the  speed 
of  a  computer  into  play,  computing  the  envelope  was  found  to  be 
cumbersome.  Instead,  using  the  area  obtained  by  integrating  the 
absolute  value  of  the  trace  amplitude  was  found  to  be  more 
efficient. 
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Mathematically,  the  definition  of  the  area  under  the  vertical 
component  of  the  Rayleigh  wave  that  is  now  being  used  is 


ARZ  =  /  | f  (t)  |dt 


where  f  (t)  is  the  seismogram  trace  and  the  group  velocity- 

window  we  are  examining. 

ERZ:  By  definition,  ERZ  is  proportional  to  the  total 

energy  in  the  Rayleigh  wave: 

T2 

ERZ  =  /  f 2 (t) dt . 

T 

1 


The  motivation  for  this  discriminant  can  be  found  in  the  need 
to  develop  some  easily  calculable  quantity  which  has  a  physical 
interpretation,  and  has  better  discrimination  ability  than  ARZ. 

REF:  This  quantity  is  defined  as  the  Rayleigh  energy  flux, 

the  energy  transported  across  a  unit  width  of  the  wave  guide  by 
the  Rayleigh  wave.  In  a  perfectly  elastic  medium  the  energy  flux 
out  of  any  closed  surface  containing  the  source  is  equal  to  the 
energy  injected  by  the  source  into  the  medium.  Therefore,  REF, 
when  corrected  for  geometrical  spreading,  shows  effects  of  only 
the  crustal  filtering,  i.e.,  scattering  and  absorption  losses, 
and  the  source  radiation  pattern. 

00  00 

REF  =  /  U(w)u ?!  p  (z) 

-00  0 

where 

U:  group  velocity 
Q:  normalized  radial  displacement 
W:  normalized  vertical  displacement 
W(0,u)):  observed  surface  vertical  displacement 
p:  density 


Q2 (z,u) 


+  W2  (z,u>) 


w  (0,a))dsda) 
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REF  -  W  (10)  ;  This  quantity  is  the  structure  correction  for 
the  REF,  which  takes  into  account  the  effect  of  structure  for  a 
particular  recording  station.  If  the  structure  under  the  re¬ 
cording  station  is  known,  then  the  normalized  spectral  energy 
density,  E  ,  may  be  calculated  using  Harkrider's  dispersion 
program. 

En  =  ^  /  P<z)  £q2(z,uj)  +  W2(z,oj)J  dz 

The  structure  correction  factor  can  be  defined  as 

K(u>)  =  U(a>)En(a>) 

Therefore,  we  can  write  the  Rayleigh  energy  flux  as 

REF  =  /  K(oj)  Wq  Coj )  d(x), 

"l 

where  the  surface  wave  packet  is  band-limited  (u^cuXu^)  as  a 
consequence  of  the  source  excitation  and  the  crustal  and  in¬ 
strumental  filtering. 

Program  Description 

A  flow  diagram  of  the  AR  program  is  shown  in  Figure  1.  It 
shows  what  steps  are  required  to  calculate  ARZ,  ERZ,  and  REF;  in 
particular,  what  spectra  are  calculated  and  where  the  instru¬ 
ment  response  and  noise  corrections  are  applied.  The  following 
paragraphs  give  a  modular  description  of  the  program. 

Tape  Input:  The  entire  trace  of  the  first  seismogram  is 
read.  The  program  will  handle  up  to  16000  points. 

Card  Input;  Besides  providing  the  event  identification 
information,  options  and  controls  for  the  tape  input  are  also 
given  (i.e.,  the  record  start  time  can  be  corrected). 

Preprocessed  Time  Series:  In  this  module,  the  trace  is 
prepared  for  computation  by  correcting  for  static  magnification, 
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decimating,  detrending,  tapering,  and  bandpass  filtering  from 

0 „ 02  to  0.10  cps. 

Computation: 

ARZ  -  The  area  under  both  the  signal  and  the  noise  sample 
are  calculated  using  the  QUAPR  subroutine  (Simpson's  rule  inte¬ 
gration)  ,  and  then  the  area  of  the  noise  is  subtracted  from  that 

of  the  signal. 

ERZ  —  First,  the  spectral  amplitude  and  the  phase  of 
the  signal  and  noise  samples  are  computed  using  the  subroutine 
GRTZSPEC,  and  the  phase  is  unwound.  Then  the  spectra  are  cor¬ 
rected  for  system  response.  Finally,  the  spectra  of  the  noise 
is  subtracted  from  that  of  the  signal  and  ERZ  is  calculated 
(see  definition) . 

Hjrp  —  using  the  spectra  calculated  for  ERZ,  the  Rayleigh 
energy  flux  is  computed  (see  definition) .  This  is  modified  by 
the  station  structure  corrections  as  given  by  Harkrider  s  pro 

gram. 

Output: 

Printout:  Besides  the  event  identification  information, 

several  plots  are  furnished: 

1.  Signal  sample  (time  series). 

2.  Noise  time  series. 

3.  Signal  amplitude  spectra  corrected  for  system 
response  and  static  magnification. 

4.  Noise  amplitude  spectra  corrected  for  system 
response  and  static  magnification. 

5.  Signal  amplitude  spectra  corrected  for  noise. 

6.  Particle  velocity  spectra  corrected  for  system 
and  noise  response. 

Also,  values  are  given  for  the  AR  of  the  signal,  the  AR 
of  the  noise,  the  AR  of  the  signal  corrected  for  noise  (ARZ), 
the  ERZ  of  the  signal  corrected  for  noise,  the  REF  of  the  signal 
corrected  for  structure,  and  the  Rayleigh  wave  magnitude  (My)  of 

the  signal. 

Binary  Save  Tape: 


With  the  event  identification,  the 


amplitude  and  phase  spectrum  will  be  saved  of  the  resultant 
(signal  less  noise)  and  noise. 

DATA  ANALYSIS  AND  RESULTS 

The  ability  of  AR  to  discriminate  between  explosions  and 
shallow  earthquakes  has  been  demonstrated  by  several  authors. 

The  purpose  of  this  work  is  to  improve  the  discrimination  ability 
of  AR  and  other  parameters  with  a  minimum  of  computing  effort. 

The  bulk  of  the  time  by  other  authors  and  in  this  investiga¬ 
tion  was  spent  attempting  to  apply  the  proper  corrections  for 
distance#  azimuth#  and  attenuation  in  order  to  eliminate  path 
effects  and  to  try  to  make  the  discriminant  a  diagnostic  of  the 
source  only.  This  can  be  graphically  demonstrated  if  the  scatter 
can  be  reduced  for  a  particular  event  in  a  plot  of  ARZ  versus 
distance. 

A  correction  which  has  been  found  to  be  particularly  fruit¬ 
ful  is  the  normalization  of  ARZ  at  all  stations  to  a  fixed  sur¬ 
face  wave  magnitude.  By  using  the  surface  wave  magnitude#  we  only 
need  to  deal  with  the  vertical  component  of  Rayleigh  wave  and  do 
not  need  to  examine  the  short-period  body  wave  train.  The  sur¬ 
face  wave  magnitude  is  defined  by 

Mv  =  Log  A/  +  1.66  Log  A°  -  0.18 


where 

1)  A/t  =  maximum  particle  velocity  peak  to  peak#  between 
T  =  17  and  30  seconds 

2)  A  =  distance  in  degrees. 

The  normalization  is  accomplished  by: 

(M  -M  )  =  Log (A/  )  -  Log  (A/T).  +  1.66 (LogA  -  Log  A.) 

V  V  •  1  1  1 

therefore# 


R  =  anti-log  (Mv~Mv  ) 

i 


where  Mv=  average  Rayleigh  magnitude 
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for  the  set  of  records 


analyzed,  i.e.,  M  =  1/n 


n 

l 

i=l 


M 


Therefore,  R  •  AR  =  normalized  value  of  AR  to  a  given  Mv* 

Figures  2  and  3  give  illustrations  on  the  application  of  this 
normalization  to  an  earthquake  (Fallon)  and  an  explosion  ^Shoal) • 
The  uncorrected  plot  of  ARZ  versus  distance  for  10  stations  is 
shown  in  Figure  2a  with  the  normalized  values  of  ARZ  versus  dis¬ 
tance  shown  in  2b.  The  reduction  in  the  scatter  of  points  due  to 
the  normalization  is  quite  obvious. 


In  an  effort  to  place  a  physical  significance  on  the  measure¬ 
ments,  the  parameter  ERZ  was  developed  (see  definition).  The  un¬ 
corrected  plot  of  ERZ  versus  distance  for  the  Fallon  earthquake 
is  shown  in  Figure  2c,  and  the  normalized  plot  is  shown  in 
Figure  2d.  Again  the  reduction  in  scatter  of  the  points  is  quite 
obvious.  The  normalization  of  ERZ  values  is  accomplished  by  the 
following: 

since  A2  ~  E  (i.e.,  Amplitude2  ~  Energy) 


and 


Therefore, 


Log (A/T)  -  Log (A/T) +3.32  (LogA-LogAi) 


anti-log  2 


M  -M 


and  ERZ  •  R_  =  normalized  value  of  ERZ  to  an  M  . 

E  v 

A  similar  set  of  figures  is  shown  for  the  Shoal  explosion. 

The  uncorrected  and  normalized  values  for  ARZ  and  ERZ  is  shown 

in  Figure  3a  through  3d.  The  reduction  of  point  scatter  is 

again  quite  significant. 


With  the  large  reduction  in  point  scatter  possible  due  to 
the  normalization  to  M  ,  further  reduction  in  the  scatter  was 
thought  possible  with  the  inclusion  of  a  station  structure 
correction.  The  parameter  REF,  the  energy  transported  across 
a  unit  width  of  the  wave  guide  by  the  Rayleigh  wave,  was 


-6- 


developed  at  this  point  (see  definition) .  Figures  4  and  5 
illustrate  the  use  of  REF  and  the  effect  of  the  structure 
correction.  Figures  4a  through  4f  deal  with  a  Guerrero, 

Mexico,  earthquake.  (For  identification  information  on  each 
event  discussed  see  Table  I  -  Event  Epicenter  Listing) .  Plots 
similar  to  those  for  Fallon  and  Shoal  of  ARZ  and  ERZ  normalized 
to  M  -  4.0  are  given  in  Figures  4a  and  4b.  In  Figure  4c,  the 
REF  normalized  to  My  =  4.0  is  plotted  versus  distance  with  no 
structure  correction  for  seven  stations.  The  standard  deviation 
of  a  least  squares  straight  line  is  SD  ='  +  0.077.  In  Figure  4d 
we  have  the  sample  plot  except  that  a  structure  correction  has 
been  applied  to  the  same  seven  stations.  The  standard  deviation 
of  a  least  squares  straight  line  is  SD  =  +  0.068.  A  comparison 
of  the  results  from  Figures  4c  and  4d  imply  that  the  structure 
correction  is  of  secondary  importance,  since  the  overall  re¬ 
duction  in  point  scatter  is  quite  negligible.  If  one  is  willing 
to  say  it  was  possible  that  the  structures  for  the  stations  used 
were  similar,  then  only  an  examination  of  Figures  4e  and  4f  is 
needed.  For  stations  FK-CO  and  UBO,  plots  of  energy  density 
versus  period  are  shown  together  with  group  and  phase  velocity 
curves.  The  energy  densities  are  quite  dissimilar.  Therefore, 
one  can  only  conclude  that  large  perturbations  in  structure 
have  a  negligible,  if  noticeable,  affect  on  the  REF  parameter. 

A  similar  set  of  plots  for  a  Komandorsky  Islands  earthquake 
of  8  February  1965  shows  the  same  results.  Figures  5a  and  5b 
give  normalized  plots  of  ARZ  and  ERZ.  The  normalized  REF  with  no 
structure  correction  is  shown  in  Figure  5c,  with  the  standard 
deviation  for  a  least  squares  fit  SD  =  +  0.094.  The  normalized 
REF  with  structure  correction  is  shown  in  Figure  5d,  with  a 
standard  deviation  SD  =+  0.115.  For  this  event,  the  structure 
correction  produced  worse  results  than  the  REF  computation  with 
no  structure  correction.  This  result  further  demonstrates  the 
lack  of  effectiveness  present  in  the  structure  correction. 

Figures  6  through  11  give  results  similar  to  those  discussed 
for  Figures  2  through  5.  They  cover  the  results  for  three  shots 
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(Bourbon,  Long  Shot,  and  Bilby),  and  for  three  earthquakes 
(Andreanof  Islands,  Calif ornia-Nevada  Border,  Komandorsky 
Islands  of  14  February  1965). 

Composite  plots  of  the  data  for  the  events  previously  dis¬ 
cussed  are  given  in  Figures  12  through  15-  Figures  12  and  14 
show  ARZ  and  ERZ  respectively,  taken  at  a  standard  distance  of 
1000  km  and  plotted  versus  body-wave  magnitude  (rr^) .  The  separa¬ 
tion  between  explosions  and  shallow  earthquakes  is  what  we 
expect  from  previously  published  results,  with  distinct  groupings 
of  the  data.  But  the  results  shown  in  Figures  13  and  15  do  not 
have  this  distinct  separation.  In  these  figures,  ARZ  and  ERZ, 
taken  at  a  standard  of  1000  km,  are  plotted  versus  surface  wave 
magnitude  (M  ) .  Although  the  explosions  generally  have  a  lower 
ARZ  or  ERZ  for  a  particular  magnitude,  there  is  no  separation  as 

before. 

An  explanation  for  this  lack  of  separation  can  be  found  if 
we  examine  our  definition  of  Mv  and  look  in  closer  detail  at  the 
spectra  for  each  event.  As  defined  previously,  My  is  calculated 
by  picking  the  maximum  amplitude,  peak  to  peak,  which  occurs 
between  T  =  17  and  30  seconds.  Hence,  we  are  averaging  the  long 
period  energy  over  the  entire  interval.  By  doing  this,  we  have 
lost  some  information  generally  inherent  in  the  wave  trains  for 
explosions  and  earthquakes.  Because  of  the  lack  of  shear  wave 
radiation  produced  by  explosive  sources,  the  energy  in  the  long 
period  Rayleigh  waves  should  be  reduced,  with  most  of  the  energy 
contained  in  the  shorter  periods.  For  earthquakes,  however, 
more  energy  should  be  contained  in  the  longer-period  Rayleigh 

waves. 

Finally,  an  attempt  has  been  made  to  extend  the  use  of 
ARZ  and  ERZ  to  low  signal-to-noise  ratios  with  the  use  of  a 
matched  filter.  The  ratio  of  ARZ's  or  the  square  root  of  the 
ratio  of  ERZ's  for  different  events  occurring  at  the  same  loca¬ 
tion  corresponds  to  the  value  of  the  parameter  3  obtained  from 
the  matched  filter  analysis  (Alexander  and  Rabenstine,  1967), 
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where 


&  =  Jx(t  +t)v (t)/Jy2 (t) , 
t  t 

x(t)  being  the  test  case  and  y (t)  the  known  waveform. 

A  comparison  of  the  ratios  with  &  was  made  for  the  two 
Greenland  Sea  Events  of  18  November  1966  studied  by  Alexander 
and  Rabenstine.  These  events  occurred  approximately  40  minutes 
apart  at  the  same  location,  the  first  having  M  of  3.6  and  the 
second  having  of  4.8.  A  summary  of  the  ratios  obtained  is 
given  in  Table  II. 

Good  agreement  exists  between  the  square  root  of  the  ratio 
of  ERZ's  and  For  those  values  for  which  a  substantial  dif¬ 
ference  exists,  an  explanation  can  be  found  in  the  method  used 
to  correct  for  noise,  which  is  not  the  same  for  both.  The  ERZ 
is  corrected  for  noise  by  a  simple  subtraction  of  spectra,  while 
the  value  for  the  noise  correction  used  in  calculating  £  is  the 
rms  n (t)  .  The  ratio  of  the  ARZ's  of  the  signal  alone  (ARZS) 
also  compares  favorably  with  &,  but  the  ratio  of  ARZ's  corrected 
for  noise  are  in  poor  agreement.  The  latter  ratio  is  generally 
much  lower  in  a  station-by-station  comparison,  indicating  that 
the  noise  sample  used  to  correct  ARZ  is  overestimated  especially 
at  low  S/N  ratios,  with  the  good  possibility  of  producing  a 
negative  ARZ. 

A  graphic  demonstration  of  the  compatability  of  these  ratios 
with  £  is  shown  in  Figure  16.  A  synthetic  case  was  produced 
where  an  actual  signal  is  buried  in  a  noise  sample  to  produce 
known  signal-to-noise  ratios.  This  signal  is  then  used  as  a 
matched  filter.  These  same  signal-noise  combinations  were  then 
processed  by  the  AR  program.  This  case  also  demonstrates  that 
the  best  overall  agreement  is  obtained  using  the  square  root  of 
the  energy  ratio.  But  more  work  is  needed  to  completely  under¬ 
stand  the  effect  of  the  noise  corrections  on  the  ratios, 
especially  for  low  S/N. 


CONCLUSIONS 


On  the  basis  of  this  investigation  of  discrimination  tech 
niques  using  the  vertical  component  of  the  Rayleigh  wave  we 
conclude  that: 


1.  In  plotting  ARZ,  ERZ,  or  REF  versus  Rayleigh  wave 
magnitude  for  a  particular  event,  the  point  scatter  is 
significantly  reduced  by  the  normalization  to  a  particular  Mv; 

2.  The  structure  correction  applied  to  the  REF  calcula¬ 
tion  has  little  if  any  effect; 


3.  Plots  of  ARZ  and  ERZ  versus  body  wave  magnitude  yield 
greater  discrimination  between  explosions  and  shallow  earth¬ 
quakes  than  plots  using  the  surface  wave  magnitude; 


4.  The  matched  filter  can  be  used  to  extend  the  ERZ 
parameter  to  low  signal-to-noise  ratios  with  reasonable  agreement 
between  the  two  parameters  (/ERZ2yERZ  and  a)  where  they  can  be 
compared. 


RECOMMENDATIONS : 

The  parameter  which  yields  the  greatest  discrimination  is 
ERZ.  In  plots  versus  body-wave  magnitude,  the  separation  is  of 
the  same  order  as  those  using  ARZ.  Although  it  takes  slightly 
more  time  to  calculate  than  ARZ,  this  is  more  than  offset  by  the 
greater  accuracy  in  the  noise  correction,  the  good  agreement  in 
comparisons  with  the  matched  filter,  and  the  fact  that  this 
parameter  his  physical  significance.  Although  the  REF  parameter 
had  a  more  satisfying  physical  interpretation,  the  negligible 
improvement  when  the  structure  correction  is  applied  leaves 
something  to  be  desired  for  the  amount  of  analysis  and  computer 
time  needed.  Therefore,  ERZ  should  be  used  as  a  basis  for 
discrimination  programs  and  further  investigations  into  the 
discrimination  problem.  These  investigations  should  include 
an  examination  of  the  matched  filter  and  the  ERZ  and  ARZ  ratios 
for  low  S/N. 
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TABLE  II 


Comparative  Data  From  Matched  Filter  and  AR  Programs 

Greenland  Sea  Events  18  Nov.  1966 


Station 

Code 

ARZ. 

ARZ2 

ARZS. 

ARZSZJ 

HN-ME 

0.0074 

0.2380 

RK-ON 

0.0445 

0.0663 

PG-BC 

0.0528 

0.0851 

KC-MO 

0.0349 

0.0711 

MO- ID 

0.0458 

0.0865 

UBO 

-0.0700 

0.0700 

AX2AL 

0.0234 

0.0805 

EU2AL 

-0.0006 

0.0870 

V'JMO 

0.0443 

0.0728 

BE-FL 

0.0321 

0.1765 

JE-LA 

0.0025 

0.1021 

MN-NV 

0.0459 

0.1586 

Notes:  Subscript  1: 

Greenland  Sea 

Subscript  2; 

Greenland  Sea 

ARZi/ARZ2: 

Ratio  of  ARZ 1  i 

ARZS1/ARZS2: 

Ratio  of  ARZ1; 

ERZ1/ERZ2: 

Ratio  of  ERZ1; 

Matched  Filter 

A 

a 

ERZ. 

ERZ  2 

0.0370 

0.0813 

0.0655 

0.0557 

0.0775 

0.0708 

0.0535 

0.0458 

0.0720 

0.0656 

0.0330 

0.0  '12 

0.0465 

0.0600 

0.0490 

0.0448 

0.0565 

0.0510 

0.0370 

0.CH75 

0.0465 

0.0837 

0.0430 

0.0960 

Eq.  18:07:54.0  Z 

Eq .  18:48:43.9  A 
corrected  for  noise 

not  corrected  for  noise 

corrected  for  noise 

EARTHQUAKE  GUERRERO  MEXICO 

13  APRIL  1967,  m  =  5.6,  h  =  86  km,  MN-NV ,  A  =  2811  km 
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Figure  1.  Flow  chart  for  Rayleigh  wave  discrimination  program 
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Figure  13.  ARZ  (At  1000  km)  vs.  Rayleigh  wave 
magnitude  (M^)  for  13  events. 
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Figure  16.  Comparison  of  data  for  synthetic  case  using 
matched  filter  and  AR  program  -  measured 
signal  level  vs.  signal  input  level  (S/N) . 
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II  SUPPLEMENTARY  NOTt* 


It  AMTRACT 

Discrimination  between  underground  nuclear  explosions  and 
shallow  earthquakes  using  the  vertical  component  of  the  Rayleigh 
wave  is  achieved  using  several  techniques.  This  analysis  includes 
the  previous  developed  area  under  the  Rayleigh  wave  (ARZ) ,  the 
newly  applied  total  energy  (ERZ)  and  the  total  energy  transported 
across  a  unit  width  of  the  waveguide  by  the  Rayleigh  wave  (REF) . 
Results  for  several  explosions  and  earthquakes  of  varying  magnitudes 
are  presented.  Evaluation  of  these  techniques  and  their  applicabil¬ 
ity  in  an  automated  discrimination  program  is  discussed.  An  attempt 
is  made  to  incorporate  a  matched  filter  approach  for  weak  signals. £ 
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